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Abstract 


A Mathematical technique is presented for 
Im pr o v ed analysis of a vide class of dynamic and 
aeroelastic systems characterized by several 
degree s-of- freedom. The technique enables greater 
utilization of the usual eigensolution obtained 
from the system dynamic equations by systematizing 
the identification of destabilizing and/or 
stiffening forces. Included, as illustrative 
examples of the use of the technique, are analyses 
of a helicopter rotor blade for bending-torsion 
divergence and flutter and for pitch- lag/ flap 
instability. 


Notation 


[A], [B], Inertia, damping and stiffness matrices, respec- 
[C] tively, Eq. (1) 




Pitch-lag coupling for k'th edgewise mode 
(=49/ Aq^) 

Pitch-flap coupling for n'th flatwise node 
(=A0/ 

Section lift curve slope, /rad 


^ij * Elements of the [A], [B] and [C] matrices 

C^ Viscous equivalent structural damping of k'th 

edgewise mode 

c Blade chord, in. 


Ely, El* 


Flatwise and edgewise bending stiffness, respec- 
tively, lb-in . 2 


{F(t)} Dynamic excitation force vector, Eq. (l) 


Resultant driving force for n'th degree of freedom, 
Eq. (5) 


[G(Xi)] Dynamic matrix for i’th eigenvalue, Eq. (3) 

p 

GJ Torsional stiffness, lb- in. 

Kj. Root feathering spring, in. -lb /rad. 

k. Polar radius of gyration of spar about its center, 

in. 


Vio’ 

**10 

m 


“o 


Section thickness-wise and chordwise mass radii of 
gyration, respectively, about spar center, in. 

? 2 

Section mass distribution, lb- sec /in. 

Reference mass distribution, (= 0.000776 lb-sec 2 /in. 2 ) 


*Senior Research Engineer, Rotary Wing Technology 
Group. 


[P A j] , "Stability" Force Phasing Matrices for i'th eigen- 
[%,], value, Bqs. ( 6 ) through ( 8 ) 

[^] 

, "Btiffh**#" Tore* Kmaing Matrices for i'th 
eigenvalue, Bqs. (10) thresh ( 13 ) 


k'th edgewise modal response • 
i’th flatwise modal response variable 




j'th torsional modal response variable, (j » 1 , 
for rigid feathering) 


R Rotor radius, in. 

r Blade spanwlse location, in. 

T Tension at r, lb 

t Time, sec 

jx} Vector of degree* of freedom 

yj_o . Chordwise positions forward of spar center of mass 
y 10 /h » center, quarter chord, and three-quarter chord, 

y i°i A **• 

a 0 Spanwiae variable section angle of attack about 

which perturbations occur, rad 


/} Blade pro-coning angle, rad 

Angle defined in Fig. 1 (» arg Aj ) 

7^ k'th assumed edgewise mode shape 

y wi i'th as maned flatwise mode shape 

y e , j'th assumed torsion mode shape 

S Qt 5^ Coefficients describing quartic variation of profile 
drag coefficient with angle of attack 

Kronecker delta 

77 Humber defined in Eq. ( 9 ) 

8 Geometric (collective) pitch angle at r, rad 

0 e Elastic torsion deflection at r, rad 

A (Uhiform) rotor Inflow 

Aj i'th eigenvalue, /sec 

O f Ji 

p Air density, lb-sec /in. 


cr 


°i 

a 


Blade solidity 

Real part of i'th eigenvalue, /sec 
i'th eigenvector of dynamic matrix equation 
Rotor rotational speed, rad/sec 
Imaginary part of i'th eigenvalue, /sec 
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(*) Differentiation with respect to (ftt) 

( )’ Differentiation with respect to Y 

(“) Indicates quantity is nondimensionalized using 

combinations of R, diq and ft, as appropriate 

£ 3 Diagonal matrix 


I. Introduction 

Dynamic and aeroelastic analyses of aerospace 
structures typically involve deriving and solving 
sets of linear differential equations of motion 
generally written in matrix fonn: 

[a]|x| + [B]|x| + [C ]|x | = |F(t)| (1) 

In general, the A, B and C matrices are square 
and real- valued. A recognized hallmark of rotary 
wing and turbcmachinery dynamics is an abundance 
of nonconservative forces (usually involving rotor 
rotation speed). Consequently, the resulting 
analyses produce matrix equations of motion of 
the above type which are highly nonsymmetrical, 
and often of large orders. 

Although a large part of the dynamic analyst’s 
job involves the calculation of dynamic loads and 
stresses due to explicit excitations, the scope 
of this paper will be limited to the equally 
important eigenproblem (F(t ) = 0): 

{*} (si 


[m»| * mv mV) . CG<x t >:y 

(3) 

The eigenvalues A (= cr±ia>), which give stability 
and natural frequency information are obtained 
from the familiar characteristic determinant: 


|[a]x+[b]\+[c]| ■ o (4) 

by any of various well-established methods 
( 2 ), '3'. The "flutter” mode shapes, are 

obtained from Eq. ( 3 ) once the eigenvalues are 
known. 


This paper presents an easily implemented 
technique for the improved analysis of dynamic 
systems of the type described above . The technique 
requires a reliable eigensolution and involves 
manipulations of the given dynamic equations, 
their eigenvalues and eigenvectors. Specifically, 
the technique systematizes the identification of 
destabilizing and/or stiffening forces by the 
calculation of "force phasing matrices". Applica- 
tions of the technique to analyses of bending- 
torsion divergence and flutter and of pitch- lag/ 
flap instability of a helicopter rotor blade are 
presented. Furthermore, this paper essentially 
represents an expansion of a portion of an earlier 
paper' 4 '. 

II , Mathematical Development 

The principal function of the force-phasing 
matrix technique is to identify those force terms 
in the equations of motion which, for an unstable 
mode, are so phased by the mode shape as to be 
drivers of the motion. The technique is perhaps 
nothing more than a formalization of the intuitive 
use an experienced dynamicist would make of the 
eigenvector information. The basis of the tech- 
nique can be seen by writing any single equation 
of the set represented by Eq. (3) as the sum of 
the mass, damper and spring forces of the diagonal 
degree-of- freedom and the remaining forces acting 
as a combined exciting force. 

a + b X. + c 

nn i Y n nn i ui nn T n 

+ ^ (a rvJ X i +b n3 Xl + ° n d )tf ?i i> = 0 ( 5 ) 

v X ' 

For the usual case a^, b^ and c nn are all 
positive numbers; that is, each mass when un- 
coupled from the others is a stable spring-mas s- 
damper system. Since the root, is generally 
complex, Eq. ( 5 ) can then be interp retted as the 
sum of four complex quantities or vectors in the 
complex plane which must, furthermore, be in 
equilibrium. Assuming that the root with 
positive imaginary part is used throughout, the 
argument of the root, y± 9 is the angle by which 
the inertia force vector is rotated relative to 
the damper force vector and the damper force 
vector is rotated relative to the spring force 
vector. For an unstable root this angle will 
be less than 90 degrees. If a purely imaginary 
value is assigned to the spring force vector, 
unstable motion is assumed and it is recalled 
that the four vectors are in equilibrium, then 
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Jm, 


the real parts of the damper and inertia force 
vectors -will he negative and the driving force 
must always have a positive real part. Figure 1, 
which demonstrates this argument, shows the four 
force vectors in the complex plane for an 
unstable oscillatory mode (Re^) = o^X)) and 
for unit imaginary displacement: 




4 c nn <t>P (SPRING FORCE) 



Figure 1. Force-Vector Diagram for n’th Degree- 
Of-Freedom, i’th Mode (Oscillatory 
Instability) 

A secondary function of the technique is to 
identify those terms in the equations which, for 
any coupled mode, act as stiffness so as to 
increase the coupled frequency of the mode. 
Reference to Figure 1 shows that driving forces 
with positive imaginary parts will tend to rein- 
force the diagonal spring term and, hence, raise 
the frequency of the coupled mode. An interesting 
observation that can be made from Figure 1 is that, 
for unstable motion, the diagonal damper force 
also has a positive imaginary part. Hence, it 
tends to stiffen the (unstable) coupled mode in 
contrast to the frequency lowering effect of 
damping for stable motion. 



Figure 2. Force-Vector Diagram for n’th Degree- of - 
Freedom, Divergence Instability. 

Again, the driving force is always a positive 
real number. Furthermore, for divergences, stif- 
fening forces are by definition stabilizing; hence 
those components of the driving force which are 
negative are also those that stiffen the coupled 
mode. 

These inte rpre tat i ons of unstable motion can 
be quantitatively implemented first, by multiplying 
each of the dynamic equations (i.e. , each row of 
the equation (l)) by a quantity which make s the 
diagonal stiffness force (stiffness matrix el ement 
x displacement) become pure im a gina ry and second, 
by representing the modal vector as a diagonal 
(square) matrix. This latter operation has the 
effect of evaluating the magnitudes of the com- 
ponent dynamic equation forces without numerically 
adding them together. The resulting "stability" 
force phasing matrices are then readily written as 

[P Ai ] - n/rp cl, 3x|LA]r (6) 

[P Bi ] » (ReCn/cp^JXitBtf qJ 1 ’] (7) 

[p„ ] - <ReC (8) 

c i 

where 


Figure 2 shows the same forces as vectors 
for an unstable aperiodic mode (divergence) for 
negative unit real displacement: 


| 1 , for oscillatory instabilities 
(-1. , for divergences 


(9] 


and where the eigenvalue in the upper half plane 
is used. 


In all cases, the real parts of the above 
indicated matrix expressions give instability 
driving force information. Forces defined by 
elements of the A, B and C dynamic equation 
matrices which are phased by the mode shape so 
as to be drivers of the motion then cause the 
corresponding elements of the Pa^, and P£j_ 

"stability” force phasing matrices, respectively, 
to be positive and proportional to their strength 
as drivers. 

Stiffening driving force information is 
obtained differently for oscillatory motion and 
for aperiodic motion. Those elements of the 
dynamic equation matrices which are phased so as 
to be stiffeners of the coupled mode will cause 
the corresponding elements of the matrix expres- 
sions to be either positive imaginary for oscil- 
latory motion, or negative real for aperiodic 
motion. The resulting "stiffness" force phasing 
matrices are then expressed as: 



= Anf Tl/«P (i) JX?[A]C cp (i) 0 

(10) 


= Anf <P (i ^ 

(H) 


= Anf T)A> ( %C]F V (1) Q 

(12) 


for oscillatory motion, and: 


[l ( )i ] = £ p ( )i ] (13) 

for aperiodic motion. 

It should be stressed that these force 
phasing matrices are no more than a more system- 
atic and efficient interpretation of the all too 
often voluminous eigensolution information. The 
following sections illustrate the usage of the 
force phasing matrix technique in substantiating 
what is generally known of same rather fundamental, 
classical helicopter rotor blade instabilities. 


III. Description of Illustrative 
Rotor Blade Example 


For illustrative purposes, relatively simple 
linear equations of motion were formulated for a 
generalized untwisted helicopter rotor blade and 
then applied to a realistic nonarticulated rotor 
configuration. The blade is assumed to be oper- 
ating in an unstalled hover condition at some 
collective angle and with a built-in coning angle. 
Perturbative elastic flatwise, edgewise and torsion 


motions are assumed to occur about the preconed 
position. The resulting linear aeroelastic 
equations are fairly standard ^)> ( )• quasi- 
static aerodynamics (uniform inflow) is assumed 
and a normal mode description of the blade 
elasticity is employed. Thus, for the chosen 
configuration, two flatwise bending modes, one 
edgewise bending mode, and the rigid feathering 
degree -of- freedom are assumed. The resulting 
response vector, jx)- , consists of the quantities 
qwi, <lw2> qvi, and qg-^ whose detailed dynamic 
equations are given in the Appendix. The dynamic 
equations then comprise a set of four differential 
equations written as a 4 x 4 matrix equation of 
the Eq. (l) type. The aeroelastic degrees-of- 
freedom together with the general parameters are 
shown in Figure 4: 



Figure 4. Schematic of Nonarticulated Rotor 

Configuration and Aeroelastic Degrees- 
of- Freedom.. 

The basic configuration incorporates a 
counterweight over the outer 70 percent of the 
blade, pitch- flap coupling (determined from the 
geometry of the pushrod attachment and flatwise 
modal deflection) and pitch- lag coupling of 
arbitrary magnitude. The chordwise position of 
the counterweight and the magnitude of pitch- lag 
coupling are purposely varied in the following 
analysis to establish known blade instabilities 
in order to illustrate the phasing matrix 
analysis technique. Table I below summarizes 
the pertinent geometric and aeroelastic data 
for the rotor blade configuration: 


i 

I 

I 

i 
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IABI£ I - BIADE CHARACTERISTICS 


Radiu. 

210 in. 

Chord ( 0 . 1 E to tip) 

13.5 in. 

Tip Speed 

650 fps 

Pitch-Flap coupling, | | 

0 . 

0.188 

Hoof feathering spring rate 

3.55 X 10 6 in. -lb/ rad 

Blade caning 

2 deg. 

Airfoil: (HftCA 0012; Mch Jfa. *0): 


a 

6.0/ rad 

»o 

0.01 


3.30/rad 1 * 

(Uniform) mass distribution 

0.000776 lb- sec 2 /in. 2 

(Uncoupled ) blade natural frequencies : 


first flatwise mode 

1.092/rev 

second flatwise mode 

2.681/rev 

first edgewise mode 

1.390/rev 

rigid feathering 

3.820/rev 

(edgewise mode structural damping) 

0.01 

(critical damping) 



Flight condition (hovering) sea level, standard 

c o lle c tive angle, 0 10 deg. 

inflow ratio, A -.0601 


The normal flatwise and edgewise mode shapes 
used are shown in Figure 5* 


IV. Application of Analysis to 
Illustrative Example 

Basic Configuration 

For purposes of comparison, the data in 
Table I was used together with a collective angle 
of 10 deg, inflow ratio of -0.0601, and with zero 
counterweight chordwise offset (from the quarter 
chord), and zero pitch-lag coupling. This basic 
case is stable in all modes as is shown by the 
following list of resulting eigenvalues: 

X 1 2 = “0.504 ± i0.960 
X,1 = -0.411 ± 12.610 
Xj' 6 = - 0.027 ± 11.398 
X lf Q = -1.472 ± 13-506 

While all the aeroelastic modes represented by 
the various eigenvalues comprise responses in all 
the four degrees of freedom, they could be 
characterized as follows: mode 1 ( A^ 2 ) first 

flatwise bending, mode 2 ( A.^ 4 ) is second 
flatwise bending, mode 3 ( * 5 , 6 ) is first edgewise 
bending, and mode 4 ( A 7 ^ 3 ) * s rigid- feathering. 

Configuration With Rearward Chordwise 
Counterweight 

If the 70 percent outer span counterweight is 
artificially shifted aft so as to place the chord- 
wise section mass center at the 32 percent chord 
point, reference to the dynamic equations (A-l, 

A-2 and A- 3 ) then yields the following A, B and 
C matrices (given in E format): 

A (inertia) Msttrlx: 



Figure 5 . Spanwise Variation of Normal Mode 
Shapes* 


•2909*00 *,0000 -.0000 *, 1864*02 

-.0000 .2006-00 *.0000 .3462-03 

-•0000 -*0000 •2592-00 -.0000 

-.1864-02 .3462-09 -.0000 .1133-03 


B (Heaping) Ifetrix: 

.2887-00 .7757-01 -.4500-02 -.1254-01 
.7992-01 .1665-00 -.6166-02 -.3309-03 
.3099-01 ,1943-01 .7901-02 -.2020-02 
•1536-07 .6326-04 -.1058-03 .3ST0-03 


C (Stiffness) Jfetrix: 

.3382-00 -.3763-01 .4633-01 -.3067-00 

-.8232*03 .1422+01 .7464-02 -.7060-01 

.4633-01 -.2974-02 .4968-00 -*4396-01 

-.2300-02 -.5423-02 -.2919-03 ,1646-02 


The first two rows of the matrix equation 
cure the equations for the first and second 
flatwise bending modes, respectively. The third 
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row is the equation for the first edgewise mode, 
while the fourth is for rigid feathering. 
Correspondingly, elements in the first two columns 
of the matrices are terns multiplying the flatwise 
tending responses and their derivatives. Simi- 
larly, third and fourth column elements are terns 
multiplying edgewise bending and rigid feathering, 
respectively, and their derivatives. 


The eigensolution for these matrices (see 
Eq. (4)) reveals the configuration to be unstable 
in both divergence and flutter: 

X i 2 ~ O.4o8, -4.466 

X 3 ^4 = 0.300 til. 789 
= -0.0088 ± il.402 
X 7,8 = -0.578 ± i3.099 

Since the equations are in nondimens ional form, 
the units of these eigenvalues are per rotor 
revolution frequency or "p". Using Eqs. (6) 
through (9) the following "stability 1 * force 
phasing matrices are written for the unstable 
divergence mode: 


A Phasing Matrix, 


-• 4840-01 

- # 0000 

-.0000 

. 1941-03 


-. 3382-00 
. 1519-01 
-• 1929*01 
. 1424-02 


-.0000 

-. 3338-01 

-.0000 

-. 1944-05 


. 3123-02 
-• 1422+01 
. 6710-02 
. 1620-03 


-.0000 

-.0000 

-. 4229-01 

-.0000 


. 4408-04 

. 1476-02 

-. 3223-02 

. 6418-06 


-. 1113-02 

-. 3316-02 

-. 4989-00 

. 4341-05 


. 5061-03 

-. 1726-02 

-.0000 

-. 1985-04 


. 9272-02 

. 4022-02 

. 5540-01 

. 1375-03 


, 1 * 953-00 

.£ 110+01 

, 2936+01 

-. 1646-02 


B Phasing Matrix, P Bl 

-. 11784-00 -. 1715-02 
-. 6016-00 -. 6791-01 
-. 5265-00 -. 1799-01 
-. 3881-08 -. 0660-06 


C Phasing Matrix, P q 1 


The larger of the destabilizing driving 
forces, which show up as positive terms, have been 
underlined for clarity. A reasonable yard-stick 
for measuring the size of the destabilizing forces 
is to compare them to the size of the stabilizing 
element in each matrix equation row. For 
oscillatory instabilities that element would be 
the diagonal damping force; for divergences, it 
would be the diagonal stiffness force. As would 
be expected of a divergence instability, the 
major destabilizing forces are displacement de- 
pendent (i.e., appear in the C matrix). By making 
additional reference to Eqs. (A-l) and (A-2), 


(A- 3) and to their evaluation given above, the 
following interpretation can be drawn from these 
results : 


1. The unstable mechanism involves a 
coupling mainly between first flatwise bending 
and rigid feathering. The mode shape, ' = 
(0.619, 0.0336, 0.0149, and 1.0), confirms this 
result . 


2. The position of the chordwise mass 
center behind the elastic axis, as indicated by 
negative dynamic equation elements a^4 and c 4^1 
and reference to the explicit statements of the 
equations in the Appendix, is a major link in 
the unstable coupling chain of events. This 
result confirms well-known results concerning the 
divergence of rotor blades ^). specifically, 
that the torsion modes drive the flatwise modes 
aerodynamic ally (elements Cj_ 4 and 4) while the 
flatwise modes drive the torsion modes with 
centrifugal inertial forces through the rearward 
mass center position (elements 04^ and c^ 2 ). 

3. The first edgewise bending mode is being 
driven by the rigid feathering through aerodynamic 
and inertia terms (element b 3 4) but is not 
actively participating in the unstable mechanism. 


In a similar manner the following force 
phasing matrices are written for the unstable 
oscillatory (flutter) mode, X^: 


A Phasing Matrix, P^ 

-. 3125-00 -.0000 
-.0000 -. 2155-00 

-.0000 -.onoo 

. 1421-02 -. 1114-04 


B Phasing Matrix, P B ^ 


C Phasing Matrix, P c ^ 

-. 1260-00 . 6779-02 

-. 4997-02 -. 1325-00 
. 3014-01 , 1173-02 

-. 6596-03 . 2794-04 


-.0000 -. 1904-01 

-.0000 -.1101-01 

-. 2731-00 -.0000 
-.0000 -. 1217-03 


-. 1968-01 

, 1282-01 

-. 4186-05 

-. 6027-03 


-. 4396-02 . 8548-00 

. 1751-01 -. 1502-00 
-.0000 . 3956-00 

-. 3155-04 -.0000 


-. 5163-00 
. 5947-00 
-. 1493-00 
• 5201-09 


. 7587-02 

-. 2979-00 

, 1659-02 

-. 5155-05 


. 2749-02 

-. 1561-01 

-. 1413-01 

-.6970-05 


Again the major destabilizing terms have 
been underlined for clarity. With few exceptions 
the same interpretations can be made of the 
flutter force phasing matrices as were made for 
the divergence ones. While the feathering 
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degree-of- freedom again drives the first flatwise 
mode aerodynaxnically (element c^_ 4), the flatwise 
mode now drives the feathering degree-of-freedam 
with vibratory inertia forces (element . 

Again these results confirm well-known findings. 


B Haasing Matrix, P B ^_ 


-. 3922-00 , 5003-02 

. 2114+01 -. 2204*00 
, 1546-01 -. 3974-03 
-. 3235 - TT 7 . 6967-05 


-. 2424-01 . 9736-02 

, 2079-00 -. 4936-02 
-. 1327-01 -, 2946-03 
, 4427-03 -. 4461-03 


Configuration With Pitch-Lag Coupling 

Using the reconfirmed knowledge that an aft 
chordwise center of mass is destabilizing, the 
configuration is altered back to the original 
quarter chord balanced configuration. In 
addition, unit pitch- lag coupling (Ay-^ = 1. ) is 
introduced into the configuration. The resulting 
dynamic equations are as follows; 

A (inertia) Matrix 


C Phasing Matrix, 


Pc 5 


. 1260-08 

. 2726-02 

. 2015-01 

-.000O 


-.<+ 643-00 

. 5299-08 

. 7210-04 

-.0000 


. 3315-00 

- TT 5 I 5 T 51 

-.0000 

-.0000 


. 7077-01 

-. 2069-00 

-. 1373-01 

-.0000 


By referring to the explicit dynamic 
equations given in the Appendix, the following 
observations can he made : 


. 2909-00 -.0000 
-.0000 . 2006-00 
-.0000 -.0000 

-.0000 -.0000 


-.0000 

-.0000 

. 2542-00 

-.0000 


-.0000 

-.0000 

-.0000 

, 1133-03 


B (Damping) Matrix 


1. The instability appears very similar to 
classic pitch-lag instability and is m a in ly a 
three-way coupling between first flatwise and 
edgewise bending modes and the rigid feathering 
degree-of -freedom. The resulting coupled mode, 
£(3) = (-0.383 - iO.435, 0.015 + 10.024, 1., 
-0.100 - i.0317). 


. 2867-00 

. 7992-01 

. 3015-01 

. 1536-07 


. 7757-01 

, 1665-00 

. 1953-01 

. 6326-04 


-. 1622-01 

-. 3596-02 

. 1003-01 

. 3801-03 


-. 1256-01 

-. 3309-03 

-. 2169-02 

, 3370-03 


C (Stiffness) Matrix 


2. The edgewise bending mode is being driven 
by inertia forces generated by flatwise bendi n g 
motion: coriolis forces proportional to precone 

and flatwise bending rate and forces proportional 
to pitch angle and flatwise bending deflection. 


, 3382-00 

-. 8232-03 

. 4633-01 

-.0000 


-. 5763-01 

, 1422+01 

-. 2974-02 

-.0000 


-. 2563-00 

-. 5358-01 

. 4431-00 

-.0000 


-. 3049-00 

-. 7113-01 

-. 4337-01 

. 1646-02 


3. The flatwise bending mode is being 
driven by aerodynamic forces generated chiefly 
by pitch- lag coupled edgewise ben d i n g and to a 
lesser extent rigid feathering deflection. 


The eigenvalues for this configuration reveal the 
configuration to he unstable in the edgewise 
bending mode: 

4 1,2 = -0.573 i 10.977 
A 3 4 = -0.4o8 ± i2.609 
45*6 = 0.0119 ± il.324 
^ 7,8 = ± 13-505 


4. The rigid feathering degree-of-freedom is 
being driven principally by a centrifugal force 
moment involving chordwise mass radii of inertia, 
pitch angle, and edgewise bending rate. 

The stiffness force-phasing matrices for this 
mode are foimed and the significant terns for the 
edgewise bending equation are underlined: 


Again, the stability force-phasing matrices 
are formed for the unstable mode ( Ac;) and the 
larger positive terms are underlined: 


A Phasing Ifetrix, P^ 


-. 9137-02 

-.0000 

-.0000 

-.0000 


-.0000 

-. 6302-02 

-.0000 

-.0000 


-.0000 

-.0000 

-. 7965-02 

-.0000 


-.0000 

-.0000 

-.0000 

-. 3558-05 


A 

A Phasing Matrix, 


-. 5097-00 

.0000 

.0000 

.0000 


.0000 

-. 3516-00 

,0000 

.0000 


.0000 
.0000 
-.4455-00 
“fo5o0 


.0000 

.0000 

.0000 

-. 1935-03 
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B Phasing Matrix, 


^3 


,3425-02 

-.3695-00 

.1722-01 

.'1469-0? 


.7*27-03 

.1975-02 

-.6235-03 

-.1570-05 


.2799-01 

-.3394-00 

.1190-03 

-.1450-02 


.37BO-02 

-.1230-02 

-.9064-03 

.3998-05 


C Phasing Matrix, 


.3382-00 .2503-02 .2923-00 

.1647-01 ,142?*0l -, 1 1 79+01 

-.1777- 01 -.4505-04 . 1*431-00 

.0000 .0000 ,0000 


-.1596-00 

.7995-00 

.4327-02 

.1646-02 


It can he seen that the principal stiffening 
terns are, not unexpectedly, the diagonal mass and 
stiffness terms. The only other significant 
stiffening terms are those involving flatwise 


bending rate and deflection which are al so the 


drivers of the unstable edgewise motion. 


That the flatwise bending deflection term is 
negative and numerically greater than the rate 
dependent term can be appreciated by noting that 
the unstable coupled edgewise mode frequency, 
1.324, is lower than the original corresponding 
stable mode frequency, I.398. 

V. Concluding Remarks 

The "force-phasing" matrices technique 
provides yet another tool for understanding 
dynsmic/aerodynamic phenomena. While it does 
not, by itself, indicate stability levels such 
as are provided by the eigensolution, it does 
complement the eigensolution by giving insight 
into the details of the dynamic configuration 
which are not directly available from the 
eigenvalues and eigenvectors alone. Moreover, 
the technique requires, in particular, eigenvector 
information as a starting point. Hence, it is 
inherently incapable of answering the more 
fundamental question of why, for any one mode, 
the eigenvector elements are indeed phased as 
they are. It should also be stressed that the 
technique is a tool to be used with, and in 
support of, engineer/analyst judgement; the 
results have to be interpreted properly, generally 
in the context of the specific application. 
Finally, the relative simplicity of the formula- 
tion makes the incorporation of the technique in 
any aeroelastic eigensolution program a straight- 
forward and easily implemented task. 


Appendix - Details of Ifrnamlc Equations 

The linear dynamic equations used to 
represent the aeroelastics of the rotor blade in 
hover are formulated using an assumed modal 
approach; the derivation is standard and uses 
the nomenclature of Reference 6. The lineariza- 
tion and subsequent simplifications are based 
upon the following assumptions: 

1. quasi- static, incompressible, nonstalled 
airloads . 

2. coincident spar center, shear center 
and tension center. 

3. zero twist. 

4. two flatwise bending modes, one edgewise 
bending mode and the rigid feathering 
degree-of- freedom. 

5. normal uncoupled bending mode shapes 
(zero twist and pitch angle). 

The flatwise bending equations are then written 
as: 


J^WiYwfc) 'itn + (5 yiOcgYwiYe j ) 

+[2By w 1 (8YV yi0 og sin8 ^) ] V 

‘ (2 5 ^o 0080 Y 8j Y w i ) (A_1) 

+ 5 sin26 VwJ 

+(®V Wl Y Vk sin e oo S e) q^ + [Sy Wl (y 10cg Ooa29 

-rP sin6)YQj] <3^ + Y^C-r 2 ^ ^ 

+ ta + V*vk) + + yio 3c /4(Ye 3 q8j 

+ \nSn + - *oYvAJl} *r - 0 
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where : 


The edgewise equations are written as: 


% n 

-[2 Sy^PYvfc + yiOcgSlne YvjJ] <W 

-[2 5(Py 10cg y Vjt + Y^) y Qj ] ^ 

+(2Y v ^Y^sin9 cos9) q^ (A-2) 

+ [ EI z Y Tfc Y T n + rrY -7k Y Tn " m cos2s Y Vi Y vn-‘ 


a a = e + xT? 


(A-4) 
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